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SUMMARY 


Most hyperelastic materials are treated as incompressible or nearly 
incompressible in analytical approach. The use of the penalty function, to 
account for near incompressibility is discussed and compared to that of 
Lagrange multiplier. A scheme to use Lagrange multiplier, without having 
to treat it as unknown, is also presented. 


INTRODUCTION 


The finite element analysis of hyperelastic materials involves large 
displacements and nonlinear material description. These materials are 
often considered as incompressible or nearly incompressible in theoretical 
developments. The incompressibility condition leads to certain simplifications 
in the analysis for exact solutions. Such is not, however, the case in 
finite element applications. 

In the variational formulation of the finite element problems of 
incompressible media, the incompressibility condition is introduced through 
the use of Lagrange multiplier (1,2). This multiplier is an additional 
unknown scalar function which can be accommodated in the discrete model 
by displacement shape functions or by similar methods. The procedure re- 
sults in an increase in the number of unknowns. 

The incompressibility condition, aside from being inconvenient in 
finite element analysis, is an approximation for rubberlike materials. 

Such approximation becomes increasingly less accurate as the percentage 
of carbon black increases in the rubber compound (3). The exact enforcement 
of incompressibility is therefore not actually required. 

In this work the near incompressibility is accounted for by use of the 
penalty function (4) in the expression for the strain energy function. 

The constitutive equations are then obtained in incremental form. This 
is most suitable for a nonlinear finite element. The incompressibility 
is shown to be closely satisfied for large penalty numbers. The 
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constitutive relations with penalty number are compared with the exact 
solution for the incompressible case. It is shown that for a given accuracy 
of the solution, the penalty number must vary in the course of finite 
element solution from element to element and from one increment to the 
next. Some numerical results are presented to illustrate these points. 

A simple scheme to use variable penalty number is also proposed. Finally 
the incremental form of constitutive equations for the case of incompressible 
material is derived from the strain energy function with the Lagrange 
multiplier. An alternative procedure is then proposed which does not require 
treatment of the Lagrange multiplier as unknown. 


BASIC EQUATION 


The strain energy function W for isotropic hyperelastic materials 
has the following form: 


w = w di, i 2 , i 3 ) 


CD 


Where Ti* I 2 , and l3are the three invariants of Green-Lagrange strain tensor 
and are related to strain components, for the axisymmetric case, as follows: 

11 = 3 + 2 ( Yu 4 - Y22 + Y33) ( 2 ) 

1 2 = 3 + 4 (^11+ S' 22 + % 33 + Y11Y22+ Y 1 1 S 3 3 + 

( 3 ) 
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1 3 = | (2 fc Yi 1 + 1) (2^22 + 1) - 4 t Yi2 2 |(l + 2 C y 33 ) (4) 


Y 11 ? Y22and Y33 are the current values of strains. In this work, only the 
axisymmetric case is considered. The discussion, however, can be extended to 
the three-dimensional case without any difficulty. In the axisymmetric case, 
the strain invariants are related to current displacements by 
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X2 is the axis of symmetry and X1X2 is the axi symmetric plane. The 
current displacements parallel to xi andx2 axes are denoted by t ui and U2 
respectively. 

t ii 

The second Piola-Kirchof f stress tensor components T at current 
time, can be obtained from the following relations 


t T « . . 1>2>3 




(9) 


To obtain the incremental form of constitutive equations we define 


ii t+At ii t ij . . 0 0 

x J = T J — T J i,j = 1,2,3 


t+At t . . -,00 

y u ■ Y u - y « i>j ■ 1,2,3 


(10) 

( 11 ) 
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where T^an d Y - are incremental stress and strain components respectively. 
The incremental stress-strain relations may now be obtained from the 
preceding equations, after some algebraic manipulations, as 


t ■ 


i» j » • * = 1,2,3 


( 12 ) 


where 


J ijkl 


3 2 W \ 3I n 31, 


3l 31 ) 3y .. 3y. . 
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3l ‘3y, . 3y . . 

m ' kl 1 x j 


j • • • 1*2,3 


(13) 


Repeated indices imply summation convention. 

STRAIN ENERGY AND PENALTY FUNCTION 


The strain energy function for nearly incompressible materials can be ob- 
tained by a series expansion of the strain energy function about (I 3 -l) 
and retaining the leading terms 

W = W 1 (I 1 ,I 2 ) + W 2 (I 2 ,Ii) (I 3 -l) 

+ W 3 (Ii,I 2 ) (I 3 -l ) 2 + . . . (14) 

We only consider the following special case 


W = Wi(Ii,I 2 ) + Hi (I 3 -I) + H 2 (I 3 -I ) 2 (15) 


where H 1 and H 2 are constants. Many different strain energy functions 
have been proposed in the literature by further expansion of W x and 
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retaining the leading terms. A general form covering the proposed models is 
w = ^ =0 C di-3) 1 (l2-3) j + Hi (I 3 -l) +H 2 (I 3 -1) 2 (16) 


where CL. are constants. The constant Hj is "not independent if the undeformed 
state is J stress free, and should satisfy the following relation 

Hi = - (C10 + 2 C01) ( 17 ) 

Co 0 = 0 

We may now consider H 2 as a penalty number to handle the imcompressibility . 

The satisfaction of incompressibility requires H 2 to approach infinity. 

For practical purpose the incompressibility can be approximately satisfied 
by not too large values of H 2 . The incompressibility can, however, be 
satisfied more accurately as H 2 gets larger. 

In finite element analysis, however, the large values of H 2 can 
lead to computational problems due to overriding stiff ness associated with 
H 2 , as discussed in (4). A scheme to employ the variable H 2 , depending 
on the local deviations from ideal incompressibility, can therefore improve 
the solution as discussed later. On the other hand we can recover the 
classical approach by letting H 2 be zero and treating Hi as the unknown 
Lagrange multiplier 


w = z: c di-3) 1 (i 2 -3) j + X(i 3 -i) 


i 3 - 1 = 0 


(18) 


Unlike the penalty function, X is then considered as an additional unknown 
and is equivalent to hydrostatic pressure. 


ONE DIMENSIONAL STRESS-STRAIN 
RELATION 


To compare the expressions (16) and (18) and to see the behaviour of the 
material with penalty number, we consider the case of one dimensional stress- 
strain relation. Let us consider an axLsymmetric medium subject to the 
following uniform strain field: 
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Yu = Y = constant 
Y 22 = Y 

Y3 3 = Y (19) 

where Y and Y are constants. Let us further assume the simplest form of 
Wj such that 


Wi = Ci (Ix-3) + C 2 (I 2 -3) 


( 20 ) 


For strain energy function (18), the stress-strain relations are 


Y is. 


T 


1 1 


T 


22 = 


T 


3 3 


T 


1 2 = 


however 


2Ci + 
2Ci + 

2Ci + 


4C 2 (1 + Y + Y) + 2A (1 -f 2 y) (1 + 2y ) 
4C 2 (1 + 2y) + 2X (1 + 2 y) 2 

4C 2 (1 + y+Y) + 2X (1 + 2y) (1 + 2“ ) 


0 


, not independent and is related to Y by 


( 21 ) 


(1 + 2 y ) (1 + 2 y ) 2 - 1 = 0 ( 22 ) 

11 3 3 

We choose A , so that T and T are both zero to simulate the uniaxial 
loading. This choice would then lead to the following stress-strain relation 
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1 - (1 + 2y) 


' ^2Ci 


+ 2 C 2 (l + 2y ) 


2 2 2 2 

The Cauchy stress O i s then related to T by 


a 22 = (1 +2y ) T' 


The one-dimensional stress-strain relation for the case of penalty 
number is now obtained from (16) 


T 11 = 2Ci + 4C 2 (1 + y + y) + 2 (1 + 2y) (1 + 2y)z 


+ 2H 2 (I 3 -l) 


T 22 = 2Ci + 4C 2 (1+27) + 2 (1 + 2y) 2 Hi + 2H 2 (I 3 -l)] 


T 3 3 = T 11 


Following the same procedure, we arrive at the following one dimensional 
constitutive equation 


1 _ 0- + 2y) 
(1 + 2 Y ) 


2Ci + 2C 2 (1 + 2y) 


where Y is now related to H 2 by 
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2H 2 |a + 2y) 2 (1 + 2y) - 1 J = 

4Ci ( Y + Y + 2 yy) + 4C 2 (Y + y + 16 Y Y ) 

2 (1 + 2Y) (1 + 2y) 


( 27 ) 


The above equation must be satisfied for all values of H It can be 

seen that, for H 2 approaching the infinity, the above equation degenerates to 


(1 + 2Y) 2 (1 + 2 y) - 1 = 0 (28) 


which is the expression of incompressibility condition. In this case 
equation (26) would become identical to equation (23). 


^COMPRESSIBILITY AND PENALTY NUMBER 


We consider the case where incompressibility is to be satisfied within 
some prescribed accuracy £ ; that is we require 


(1 + 2Y) 2 (1 + 2Y) - 1 = £ 

For numerical illustration, consider the case where 


(29) 


= 4 


(30) 


Combining equations (27) to (30), we arrive at the following relation 


HS - 6 - j - n a +E > V2 


(31) 


where 
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X = 1 + 2y 


2H 2 
C 2 


( 32 ) 


and 


e « 1 


For small values of e , the relation (31) can be further simplified to 


He 


_ 1_ x + 4 

X xV 2 


(33) 


The stress-strain relation of the equation (23) is plotted in Figure (1) • 
The variation of He as function of strain is plotted in Figure (2). 

It can be seen that as Y increases, H must also increase accordingly, 
to maintain the same accuracy £ on incompressibility condition. It can 
also be noted that higher values of H are required in compression than in 
tension. The relation (26) may now be written as follows: 


T 


22 


com 


1 


(1 +e) /2 

(1 + 2 y ) 3 / 2 


^4 + (1 + 2Y) 



(34) 


Comparing (34) and (23), it is observed that the stress for a nearly incom- 
pressible model is always less than that of an incompressible model, for the 
same strain. This difference, however, depends on the £ and approaches 
zero as £ approached zero, or H approaches infinity. The relative error, 
however, is 


e 


22 _ 22 
j TLC — 


T 


22 

inc 


2e 

(1 + 2Y) 3 /2 - 1 


(35) 


For a fixed H, £ increases with y but the relative error in the stress 
is governed by (35) which is less sensitive to a variation in H . 
Eliminating £ between (35) and (33), we obtain 
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He 


(36) 



The above equation relates the magnitude of H to level of strain for 
a prescribed error e . The equation (33) or (36) may serve as an 
approximate method of updating H, in a problem of combined stresses, for 
improved accuracy in incompressibility or stress calculation. Further 
work is, however, required to develop a more vigorious scheme 
in the general case. 


INCREMENTAL FORM OF LAGRANGE MULTIPLIER 


We now consider an alternative approach to incompressibility problems. 
Let us consider the following strain energy function 


W = Wi(Ii, I 2 ) + A (I 3 -l) 


I 3 - 1 = 0 


(37) 


where A is the Lagrange multiplier. The second Piola-Kirchof f stresses 
can now be obtained as follow 




3W 

9 ^ 


9(I 3 - 1) 


( 1 3 - 1) 


9 t X 


ij 


(38) 


■% 

where A is the current value of A • The incremental stress-strain 
relations may be obtained by a procedure similar to that used in the 
preceding pages. This would lead to 
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X = 


t+At x - fc x 


(40) 


where A is the increment in the Lagrange multiplier between two consecutive 
steps in the incremental solution. The incremental form of incompressibility 
condition is 


3I 3 

aV. 

ij 


y. . = o 

9 i -i 


(41) 


multiplying both sides of (39) by Y^, and use of (41), leads to 


T y . . 
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The following scheme may now be considered for the solution of 
incompressible materials without the need for introducing the unknown Lagrange 
multiplier. In the process of incremental solution, the last term in the 
right hand side of equation (39) is ignored. At the end of each increment, 
t X can be updated from relation (43). There is, therefore, no need of 
treating ^ as an independent unknown. No numerical work, however, has 
been carried out with this alternative proposed scheme at this time. 
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